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STOCHASTIC 2D HYDRODYNAMICAL TYPE SYSTEMS: 
WELL POSEDNESS AND LARGE DEVIATIONS 

IGOR CHUESHOV AND ANNIE MILLET 

(— ^ , Abstract. We deal with a class of abstract nonlinear stochastic models, which covers 

(^ • many 2D hydrodynamical models including 2D Navier-Stokes equations, 2D MHD models 

^SJ ' and 2D magnetic Benard problem and also some shell models of turbulence. We first 

prove the existence and uniqueness theorem for the class considered. Our main result is 
fl\ ' a Wentzell-Freidlin type large deviation principle for small multiplicative noise which we 

^^ I prove by weak convergence method. 

(N 

'~7', 1. Introduction 

Ph ■ In recent years there has been a wide-spread interest in the study of quahtative proper- 

(~| ! ties of stochastic models which describe cooperative effects in fluids by taking into account 

macroscopic parameters such as temperature or/and magnetic field. The corresponding 
mathematical models consists in coupling the stochastic Navier-Stokes equations with 
some transport or/and Maxwell equations, which are also stochastically perturbed. 

Our goal in this paper is to suggest and develop a unified approach which makes it 
^ , possible to cover a wide class of mathematical coupled models from fluid dynamics. Due 

f^ I to well-known reasons we mainly restrict ourselves to spatially two dimensional models. 

Our unified approach is based on an abstract stochastic evolution equation in some Hilbert 
space of the form 
j^ , dtu + Au + B{u,u) + n{u) = a(t,u)W, (1.1) 

^D I where a{t, n) T4^ is a multiplicative noise white in time with spatial correlation. The hy- 

potheses which we impose on the linear operator A, the bilinear mapping B and the oper- 
ator 7^ are true in the case of 2D Navier-Stokes equation (where 7^ = 0), and also for some 
other classes of two dimensional hydrodynamical models such as magneto-hydrodynamic 

^ ' equations, the Boussinesq model for the Benard convection and 2D magnetic Benard prob- 

lem. They also cover the case of regular higher dimensional problems such as the 3D Leray 
a-model for the Navier-Stokes equation and some shell models of turbulence. See a further 
discussion in Sect. 2.1 below. 

For general abstract stochastic evolution equation in infinite dimensional spaces we refer 
to [12]. However the hypotheses in [12] do not cover our hydrodynamical type model. We 
also note the stochastic Navier-Stokes equations were studied by many authors (see, e.g., 
[5, 18, 28, 36] and the references therein). 

We first state the result on existence, uniqueness and provide a priori estimates for a 
weak (variational) solution to the abstract problem of the form (1.1) where the forcing term 
also includes a stochastic control term with a multiplicative coefficient (see Theorem 2.4). 
As a particular case, we deduce well posedness when the Brownian motion W is translated 
by a random element of its Reproducing Kernel Hilbert Space (RKHS), as well a priori 



OO 



o 



s 



2000 Mathematics Subject Classification. Primary 60H15, 60F10; Secondary 76D06, 76M35. 
Key words and phrases. Hydrodynamical models, MHD, Benard convection, shell models of turbulence, 
stochastic PDEs, large deviations. 

This research was partly supported by the research project BMF2003-01345 (A. Millet). 

1 



2 I. CHUESHOV AND A. MILLET 

bounds of the solution with constants which only depend on an a.s. bound of the RKHS 
norm of the control. In all the concrete hydro dynamical examples described above, the 
diffusion coefficient may contain a small multiple of the gradient of the solution. Thus, this 
result contains the corresponding existence and uniqueness theorems and a priori bounds 
for 2D Navier-Stokes equations (see, e.g. [28, 34]), for the Boussinesq model of the Benard 
convection (see [17], [14]), and also for the GOY shell model of turbulence (see [1] and 
[27]). Theorem 2.4 generalizes the existence result for MHD equations given in [2] to the 
case of multiplicative noise and also covers new situations such as the 2D magnetic Benard 
problem, the 3D Leray a-model and the Sabra shell model of turbulence. 

Our argument mainly follows the local monotonicity idea suggested in [28, 34]. However, 
since we deal with an abstract hydrodynamical model with a forcing term which contains 
a stochastic control under a minimal set of hypotheses, the argument requires substantial 
modifications compared to that of [34] or [27]. It relies on a two-step Gronwall lemma (see 
Lemma 4.1 below and also [14]). 

Our main result (see Theorem 3.2) is a Wentzell-Freidlin type large deviation principle 
(LDP) for stochastic equations of the form (1.1) with a := \fea as e — > 0, which describes 
the exponential rate of convergence of the solution u := vF to the deterministic solution u^. 
As in the classical case of finite-dimensional diffusions, the rate function is described by an 
energy minimization problem which involves deterministic controlled equations. The LDP 
result is that which would hold true if the solution were a continuous functional of the noise 
W . Our proof consists in transferring the LDP satisfied by the Hilbert- valued Brownian 
motion y/eW to that of a Polish-space valued measurable functional of \feW as established 
in [3]; see also [4], [13] and [15]. This is related to the Laplace principle. This approach 
has been already applied in several specific infinite dimensional situations (see, e.g, [34] 
for 2D Navier-Stokes equations, [14] for 2D Benard convection, [4] for stochastic reaction- 
diff'usion system, [25, 31] for stochastic p-Laplacian equation and some its generalizations, 
[27] for the GOY shell model of turbulence). We also refer to [6] for large deviation results 
for evolution equations in the case on non-Lipschitz coefficients. 

Our result in Theorem 3.2 comprehends a wide class of hydrodynamical systems. In 
particular, in addition to the 2D Navier-Stokes equations and the Boussinesq model men- 
tioned above, Theorem 3.2 also proves LDP for 2D MHD equations, 2D magnetic Benard 
convection, 3D Leray a-model, the Sabra shell model and dyadic model of turbulence. 
Note that unlike [34] and [27], in order to give a complete argument for the weak conver- 
gence (Proposition 3.4) and the compactness result (Proposition 3.5), we need to prove 
a time approximation result (Lemma 3.3). This requires to make stronger assumptions 
on the diffusion coefficient a, which should have some Holder time regularity, and in the 
explicit hydrodynamical models, no longer can include the gradient of the solution (see 
also [14]). 

Note that the weak convergence approach has been used recently to prove LDP for 
stochastic evolution equations which satisfy monotonicity and coercivity conditions by 
J. Ren and X. Zhang [31] and by W. Liu [25]. This class of models does not contain the 
hydrodynamical systems considered in this paper and the main PDE model for this class 
is a reaction-diffusion equation with a nonlinear monotone diffusion term perturbed by 
globally Lipschitz sub-critical nonlinearity. Let us point out one of the main differences 
which explains why we have to impose some more time regularity assumptions on the 
diffusion coefficient, in contrast with [31]. Indeed, unlike the situation considered in [31] 
we do not assume the compactness of embeddings in the corresponding Gelfand triple 
V' C H C V. Thus the elegant method in [31], which is based on some compactness 
property of the family of solutions in C([0, T], V') obtained by means of the Ascoli theorem, 
cannot be applied here, and we use a technical time discretization. Let us also point out 
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that due to the bihnear term which arises in hydrodynamical models, the control of any 
moment of the V norm for time increments of the solution is not as good as that in [31]. 
On the other hand, giving up the compactness conditions allows us to cover the important 
class of hydrodynamical models in unbounded domains. The payoff for this is some Holder 
condition in time for the diffusion coefficient (see condition (C4) below). The paper by 
W. Liu [25] does not assume the compactness embeddings in the Gelfand triple and uses 
another way to obtain the compactness of the set of solutions in some set of time continuous 
functions. This approach is based on the rather strong approximation hypothesis which 
involves some compact embeddings (see Hypotheses (A4) and (A5) in [25]). The technique 
introduced by W. Liu might be used in our framework, in order to avoid the extra time 
regularity of a, at the expense of some compact approximation condition of diffusion term. 
The corresponding proof is quite involved and we do not adapt it to keep down the size of 
the paper; we rather focus on the main technical problems raised by our models, even in 
the simple case of a time- independent diffusion coefficient. Note that even in the case of 
monotone and coercive equations, the diffusion coefficient considered in [31] (resp. [25]) 
cannot involve the gradient in order to satisfy condition (H5) (resp. (A4)). 

The paper is organized as follows. In Section 2 we describe our mathematical model with 
details and provide the corresponding motivations from the theory of (coupled) models of 
fluid dynamics. In this section we also formulate our abstract hypotheses and state the 
results about well posedness and apriori bounds of the abstract stochastic equation which 
also may contain some random control term. The proof of these technical results is given 
in the Appendix, Section 4. Note that these preliminary results are proved in a more 
general framework than what is needed to establish the large deviation principle. Indeed, 
we use them in [10] where we characterize the support of the distribution of the solution 
to the stochastic hydrodynamical equations. We formulate and prove the large deviations 
principle by the weak convergence approach in Section 3. There we use properties (such as 
a priori bounds and localized time increment estimates) of this stochastic control system 
as a preliminary step in order to apply the general LDP results from [3, 4] in our situation. 

2. Description of the model 



Let {H, |.|) denote a separable Hilbert space, A be an (unbounded) self-adjoint positive 
linear operator on H. Set V = Dom[A^). For v ^V set ||u|| = lA^tij. Let V' denote 
the dual of V (with respect to the inner product (., .) of H). Thus we have the Gelfand 
triple V G H C V . Let {u,v) denote the duality between u G V and v G V' such that 
{u, v) = (u, v) for u G V, V £ H, and let B : V x V ^ V he a continuous mapping 
(satisfying the condition (CI) given below). 

The goal of this paper is to study stochastic perturbations of the following abstract 
model in H 

dtu{t) + Au{t) + B{u{t),u{t)) + Ru{t) = f, (2.1) 

where i? is a linear bounded operator in H. We assume that the mapping B : V xV ^ V 
satisfies the following antisymmetry and bound conditions: 

Condition (CI): 

• B : V X V ^ V is a bilinear continuous mapping. 

• Form £V, i = 1,2,3, 

{B{ui,U2) , Us) = - {B{ui,U3) , U2). (2.2) 

• There exists a Banach (interpolation) space 7i possessing the properties 
(i)VcncH; 
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\v\\t^ < ao\v\ \\v\\ for any v (^ V; (2-3) 



fa) there exists a constant cq > such that 

|2 

\n 
(in) for every r/ > there exists C^ > such that 

\{B(ui,U2), U3)\ <7]\\u3f + Cr,\\ui\\l^\\u2\\n, forui £ V, i = 1,2,3. (2.4) 

Remark 2.1. (1) The relation in (2.4) obviously implies that 

|(B(ni,M2), ■"3)1 < C'lll'WsiP + (^2 II^iIIh ll'"2||^, for Ui £ V, i = 1,2,3, (2.5) 

for some positive constants Ci and C2. On the other hand, if we put in (2.5) rjC^ U3 
instead of U3, then we recover (2.4) with Crj = CiC2T]~^ Thus the requirements (2.4) 
and(2.5) are equivalent. If for u^ ^ we put now rj = 111*1117^11^211^11^311""^ in (2-4) with 
Crj = CiC2r]~^, then using (2.2) we obtain that for some constant C > 0, 

|(S(ni,n2), Ms)! < C||ui||7^ ||n2|| II^sIIt^, ior m £ V, i = 1,2,3. (2.6) 

It is also evident that (2.6) and (2.2) imply (2.4). Thus the conditions in (2.4), (2.5) and 
(2.6) are equivalent to each other. 

(2) To lighten notations for ui G V, set B{ui) := B{ui,ui); relations (2.2), (2.3) and 
(2.6) yield for every 77 > the existence of C,, > such that for ui,U2 £ V , 

\{B{ui) , U2)\ < r/ ll^if + C, |ni|2 ||m2||^. (2.7) 

Relations (2.2) and (2.7) yield 

\{B{ui)-B{u2) , Ui-U2)\ = \{B{ui-U2),U2)\ < r]\\ui-U2\\'^ + Crj\ui-U2\'^\\u2\\fi. (2.8) 

2.1. Motivation. The main motivation for the condition (CI) is that it covers a wide 
class of 2D hydrodynamical models including the following ones. An element of M^ is 
denoted u = (n^, n^). 

2.1.1. 2D Navier-Stokes equation. Let D he a bounded, open and simply connected do- 
main of M^. We consider the Navier-Stokes equation with the Dirichlet (no-slip) boundary 
conditions: 

dtu — vAu + uVu + 'Vp = f, divn = in D, u = on dD, (2.9) 

where u = (u^ {x , t) , u'^ {x , t)) is the velocity of a fluid, p(a;, t) is the pressure, z/ the kinematic 
viscosity and f{x,t) is an external density of force per volume. Let n denote the outward 
normal to dD and let 

iJ(^) = {f £ [L'^{D)] ^ div / = in D and / . n = on dD} 

he endowed with the usual L^ scalar product. Here above we set div/ = '^i=i2difi, 
Projecting on the space i7(i) of divergence free vector fields, problem (2.9) can be written 
in the form (2.1) (with i? = 0) in the space Hn\ (see e.g. [35]), where A is the Stokes 
operator defined by the bilinear form 

a('Ui,'U2) = z^V / Vul-Vu^dx, (2.10) 

with ui,U2£V = Vi= [H^iD)] ^ n H^^y The map B = Bi : Vi x Vi ^ V( is defined by 

{Bi{ui,U2) , Us) = / [ui{x)Vu2{x)]us{x)dx = 2, / u\dju\u\dx, Ui£Vi. (2.11) 
Jd ,^^^Jd 
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Using integration by parts, Schwarz's and Young's inequality, one checks that this map 
Bi satisfies the conditions of (CI) with TC = [L^(L')] n -ff(i)- The inequahty in (2.3) is 
the well-known Ladyzhenskaya inequality (see e.g. [9] or [35]). 

We can also include in (2.9) a Coriolis type force by changing / into / — Ru, where 
R{v}-,v?) = co(— «^,«^), for some constant cq. In this case we get (2.1) with i? / 0. 

The case of unbounded domains D (including Z) = M^) can be also considered in our 
abstract framework. For this we only need to shift the spectrum away from zero by 
changing A into A + Id and introducing R = —Id. 

2.1.2. 2D magneto-hydrodynamic equations. We consider magneto-hydrodynamic (MHD) 
equations for a viscous incompressible resistive fluid in a 2D domain D, which have the 
form (see, e.g., [29]): 



dtu - uiAu + nVn = -V f p + -|6r J + sbVb + /, (2.12) 

dtb-iy2Ab + uVb = bVu + g, (2.13) 

divu = 0, div6 = (2.14) 

where u = {u^(x, t),u'^{x, t)) and b = {b^{x, t), 6^(x, t)) denote velocity and magnetic fields, 
p{x, t) is a scalar pressure. We consider the following boundary conditions 

n = 0, b.n = 0, dib'^ - ^36^ = on dD (2.15) 

In equations above z^i is the kinematic viscosity, 1^2 is the magnetic diffusivity (which 
is determined from magnetic permeability and conductivity of the fluid), the positive 
parameter s is defined by the relation s = Ha?viV2-, where Ha is the so-called Hartman 
number. The given functions / = f{x,t) and g = g{x,t) represent external volume forces 
and the curl of external current applied to the fluid. We refer to [23], [16] and [33] for the 
mathematical theory for the MHD equations. 

Again, the above equations are a particular case of equation (2.1) for the following 
spaces and operators which satisfy (CI). To see this we first note that without loss of 
generality we can assume that s = 1 in (2.12) (indeed, if s 7^ 1 we can introduce a new 
magnetic field b := ^/sb and rescale the curl of the current g := y/sg). For the velocity 
part of the MHD equations, we use the same spaces Hn\ and Vi and the Stokes operator 
generated by the bilinear form defined by (2.10) with z^ = z^i. Now we denote this operator 
by^i. 

As for the magnetic part we set Ht2) = Hm and V2 = [H^{D)'j n Hf2) and define 
another Stokes operator A2 as an unbounded operator in Hf2) generated by the bilinear 
form (2.10) with z^ = z^2 when considered on the space V2. 

As in the previous case we can write (2.12)-(2.15) in the form (2.1) in the space H = 
i?(i) X F(2) with A = AixA2,R = 0. We also set F = Fi x I/2 and define B -.VxV ^V 
by the relation 

{B{zi,Z2),Z3) = {Bi{ui,U2),U3) - (Bi(6i, 62), ""s) + {Bl{ui,b2),b3) - {Bi{bi,U2),b3) 

for Zi = {ui,bi) £ V = Vi X V2, where Bi is given by (2.11). The conditions in (CI) are 
satisfied with H = ( [L^{D)Y x [L^{D)Y) n H. 

2.1.3. 2D Boussinesq model for the Benard convection. The next example is the following 
coupled system of Navier-Stokes and heat equations from the Benard convection problem 
(see e.g. [19] and the references therein). Let D = {0,1) x (0, 1) be a rectangular domain 
in the vertical plane, (61,62) the standard basis in M^ and x = (x^,a;^) an element of M^. 
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Denote hy p{x,t) the pressure field, f,g external forces, u = {u^{x,t),u'^{x,t)) the velocity 
field and 9 = 6{x,t) the temperature field satisfying the following system 

dtu + uVu - uAu + Vj9 = 6^62 + /, div n = 0, (2.16) 

dte + uVe-u"^ -kA9 = g, (2.17) 

with boundary conditions 

u = t 9 = on x^ = and x^ = 1, 

u,p, 9, u^i , 9^1 are periodic in x with period l} 

Here above v is the kinematic viscosity, k is the thermal diffusion coefficient. Let 
^(3) = |n G [L'^{D)] , divn = 0, n\2=o = ^^1x2=1 = 0, m^L.i=o = u^\x^=ij 

and if (4) = L?'{D). We also denote 

V3 = -^ n G i/(3) n \H (D)] , w|a;2=o = ^1x2=1 = 0, n is /-periodic in x > , 

Va = {9 e H^{D), 6*1^2=0 = 0\x^=i = 0, 6* is /-periodic in x^} . 

Let ^3 be the Stokes operator in H/-^\ generated by the bilinear form (2.10) considered on 
V3 and A4 be the operator in Hu-\ generated by the Dirichlet form 

a{9i,92) = K [ V9i-V92dx, ^1,^2 £^4- 
Jd 

Again, the above equations are a particular case of equation (2.1) for the following spaces 

and operators which satisfy (CI). Let H = H^^-^ x 7/(4) and V = V3 x V4. We set 

A{u, 9) = {A^u , A49), R{u, 9) = -{9e2 , u^), and define the mapping B : V x V ^ V hy 

the relation 

{B{zi,Z2),Z3) = {Bi{ui,U2),U3) + ^ / Ui 5j 6*2 6*3 dx 

for Zi = {ui,9i) £ V = V3 X V4, where Bi is given by (2.11). With these notations, 
the Boussinesq equations for {u,9) are a particular case of (2.1) with condition (CI) for 
n= {[L^{D)Y xL^{D))nH. 

2.1.4. 2D magnetic Benard problem. This is the Boussinesq model coupled with magnetic 
field (see [20]). As above let D = (0,/) x (0,1) be a rectangular domain in the vertical 
plane, (61,62) the standard basis in R^. We consider the equations 



dtu + uVu- uiAu + V {p+ -\b\^] - sbVb = 9e2 + f, div'u = 0, 

dt9 + uV9-u^-KA9 = /, 

dtb - 1/2 A6 + uVb - bVu = h, div 6 = 0, 

with boundary conditions 

u = k 9 = k b'^ = 0, d2b^ = on x^ = and x^ = 1, 

u,p, 9, b, n^i, 9^1, b^i are periodic in x with period /. 

As for the MHD case we can assume that s = 1. In this case we have (2.1) for the variable 
z = {u,9,b) with H = ii(3) x iJ(4) x -^(5), where if(3) and i/(4) are the same as in the 
previous example and ii(5) = H[3)- We also set 1/ = V3 x V4 x V5, where V3 and V4 are 



Here and below this means that cj>\^i^q — cj>\xi^i for the correspondmg function. 
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the same as above and V5 = Hi^\ n [i7^(D)] . The operator A is generated by the bihnear 
form 

a{zi,Z2) = ^1^2 ^^1 ■ '^"2 dx + K VOi ■ V6'2 dx + 1/2 V / V6{ • V6^ dx 
~[Jd Jd ^ Jd 

for Zi = {ui,9i,bi) G V. The bihnear operator B is defined by 

{B{zi,Z2),Z3) = {Bi{ui,U2),U3) - {Bi{bi,b2),U3) 

+ (Si (m, 62), 63) -(Si (61,^2), 63)+ V / 1.15^02^3^2; 

for Zi = {ui, 9i, hi) G V, where Bi is given by (2.11). We also set R{u, 6, h) = —(6*62 , u^, 0). 
It is easy to check that this model is an example of equation (2.1) with (CI), where 
7i = ( [L\D)\ ^ X L\D) X [L\D)\ ^)nH. 

2.1.5. 3D Leray a-model for Navier-Stokes equations. The theory developed in this paper 
can be also applied to some 3D models. As an example we consider 3D Leray a-model (see 
[24]; for recent development of this model we refer to [7, 8] and to the references therein). 
In a bounded 3D domain D we consider the following equations: 

dtu - uAu + vVu + \/p = f, (2.18) 

{l-aA)v = u, div'u = 0, div'i; = in D, (2.19) 

v = u = on dD. (2.20) 

where u = (u^jU^jU^) and v = {v^ , v'^ , v'^) are unknown fields, p{x,t) is the pressure. In 
the space 

H = {ue [L'^{D)] ^ : div n = in L> and n . n = on dD} 

problem (2.18)-(2.20) can be written in the form 

Ut+Au + B{GaU, u) = /, 



where A is the corresponding 3D Stokes operator (defined as in the 2D case by the form 

a{ui,U2) = i^Ej=i j 
Green operator and 



a{ui,U2) = V Y^]^-^ J^ Vu{ Vui dx on V = Hr] [H^{D)Y), G^ = {id + av-^A) ^ is the 



3 „ 

(i?(ui, 1x2) , U3) = V / u{ 9,4 uidx, UiGV = Hn [H^{D)]\ 

Note that the embedding H^/'^{D) C L^{D) implies that the inequality (2.3) holds true for 
T[ = ^L^{D)'j nH. Furthermore, Holder's inequality and the embedding H^{D) C L^{D) 
imply that for ui,U2,U3 E V, 

\{B{GaUi,U2) , ^3)1 < C\\U2\\ \GaUi\L6(^D) \u3\vi(D) < C\\u2\\ \\GaUi\\ \u3\l3(^d) 

where the last inequality comes from the fact that A^Ga is a bounded operator on H, so 
that IIGqUiII = l^aGaUil < C\ui\ < C\ui\i^'it£,y By Remark 2.1(1) this implies condition 

(CI) for Ba{ui,U2) := B{GaUl,U2). 
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2.1.6. Shell models of turbulence. Let H he a set of all sequences u = {ui,U2, ■ ■ ■) of 
complex numbers such that ^^ l^nP < oo. We consider H as a real Hilbert space endowed 
with the inner product (•, •) and the norm | • | of the form 

oo oo 

{u,v) = Re^Unvl, |np = ^|'u„P, 

n=l n=l 

where u* denotes the complex conjugate of Vn- In this space H we consider the evolution 
equation (2.1) with R = and with linear operator A and bilinear mapping B defined by 
the formulas 

f °° 1 

{Au)n = i^k'^Uu, n = 1, 2, . . . , Dom{A) = <u e H : ^J ^nknP < co > , 

where u > 0, kn = kofi"^ with feg > and fJ- > 1, and 

for n = 1, 2, . . ., where a and 6 are real numbers (here above we also assume that U-i = 
uq = v-i = vq = 0). This choice of A and B corresponds to the so-called GOY-model 
(see, e.g., [30]). If we take 

[B{u, v)]^ = -i {akn+lul^iVn+2 + bknUn-lVn+1 + «^n-l^in-l^'n-2 + bkn-lUn~2Vn^l) , 

then we obtain the Sabra shell model introduced in [26]. In both cases the equation (2.1) 
is an infinite sequence of ODEs. 

One can easily show (see [1] for the GOY model and [11] for the Sabra model) that the 
trilinear form 

oo 

{B{u,v),w) = Re y2[B{u,v)]nW^ 

n=l 

possesses the property (2.2) and also satisfies the inequality 

\{B{u,v),w)\<C\u\\A^/^v\\w\, yu,weH, \/v e Dom{A^/^). 
Thus by Remark 2.1(1) the condition (CI) holds with 7i = Dom{A^) for any choice of 

se [0,1/4]. 

We can also consider the so-called dyadic model (see, e.g., [21] and the references therein) 
which can be written as an infinite system of real ODEs of the form 

dtun + i/A2°"n„ - X^ul_, + A"+iu„u„+i = U, n = 1, 2, . . . , (2.21) 

where iy,a > 0, X > 1, uq = 0. Simple calculations show that under the condition 
a > 1/2 the system (2.21) can be written as (2.1) and that condition (CI) holds for 
[B{u,v)]n = -A"u„_iv„_i + A"+^ UnVn+1 and {Au)n = z^ A^"" 



n^ 



2.2. Stochastic model. We will consider a stochastic external random force / of the 
equation in (2.1) driven by a Wiener process W and whose intensity may depend on the 
solution u. More precisely, let Q be a linear positive operator in the Hilbert space H which 
belongs to the trace class, and hence is compact. Let Hq = Q^H. Then Hq is a Hilbert 
space with the scalar product 

together with the induced norm | • |o = vOTOo- The embedding i : Hq -^ H is Hilbert- 
Schmidt and hence compact, and moreover, i i* = Q. Let Lq = Lq{Hq,H) be the space 
of linear operators S : Hq i-^ H such that SQ^ is a Hilbert-Schmidt operator from H to 
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H. The norm in the space Ln is defined by \S\\ = tr{SQS*), where S* is the adjoint 
operator of S. The Lg-norm can be also written in the form 

oo oo 

\S\Iq = tr{[SQ'/^][SQ'/^r) = Y^ \SQ'/^i;k\^ = ^ KSQ^/YV^fcl' (2-22) 

fc=i fe=i 

for any orthonormal basis {ipk} in H. 

Let W{t) be a Wiener process defined on a filtered probability space (17,^, jFt,F), 
taking values in H and with covariance operator Q. This means that W is Gaussian, has 
independent time increments and that for s,t >0, f,g £ H, 

E{Wis),f) = and E{W{s), f){W{t),g) = {s At) {Qf,g). 

We also have the representation 

n 

W(t)= lim WJt) in L^i^-H) with W„(t) = V gy^/3,(t)e,-, (2.23) 

where (3j are standard (scalar) mutually independent Wiener processes, {ej} is an or- 
thonormal basis in H consisting of eigen-elements of Q, with Qe^ = QjCj. For details 
concerning this Wiener process we refer to [12], for instance. 

The noise intensity a : [0, T] x y — > Lq{Hq, H) of the stochastic perturbation which we 
put in (2.1) is assumed to satisfy the following growth and Lipschitz conditions: 

Condition (C2): a E C([0, T] x V; Lq[Hq, H)^ , and there exist non negative constants 
Ki and Li such that for every t G [0, T] and u,v £ V : 

(i) \a{t,u)\l^ < Ko + Ki\u\^ + K2\\u\\^ 

(ii) \(7{t,u) - a{t,v)\1 < Li\u - v\'^ + ^211^'' - vW^- 

Remark 2.2. Assume that a £ C([0, T] x Dom{A'');LQ{Ho, H)) for some s < 1/2 is such 
that 

\ait, u)\l^ <K'o + K[\A'u\\ |a(t,n) - c7it,v)\l^ < L'\A'iu - v)\^ 

for every t £ [0, T] and u,v £ Dom{A^) with some positive constants K'q, K[ and L' . By 
interpolation we have that for some constant cq > and any r] > and u £ V: 

\A'u\' < co\A'/^u\^'\u\^-^' < ^r]\A'/^u\' + Cr,\u\\ (2.24) 

Therefore in this case the conditions in (C2) are valid with positive constants K2 and L2 
which can be taken arbitrary small. This observation is important because in Theorem 2.4 
below we impose some restrictions on the values of the parameters K2 and L2. 

Recall that for u £V, B{u) = B{u,u). Consider the following stochastic equation 

du{t) + [Au{t) + B{u{t)) + Ru{t)] dt = a{t, u{t)) dW{t). (2.25) 

For technical reasons, in order to prove a large deviation principle for the law the solution 
to (2.25), we will need some precise estimates on the solution of the equation deduced 
from (2.25) by shifting the Brownian W by some random element (see e.g. [34] and [14]). 
This cannot be deduced from similar ones on u by means of a Girsanov transformation 
since the Girsanov density is not uniformly bounded when the intensity of the noise tends 
to zero (see [14]). Thus we also need to consider the corresponding shifted problem. 

To describe a set of admissible random shifts we introduce the class A as the set of 
ffo— valued (.Fj)— predictable stochastic processes h such that Jq |/i(s)|Q(is < cxd, a.s. Let 

SM = [h£ L\0,T;Ho) : / \h{s)\lds < m}. 
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The set Sm endowed with the fohowing weak topology is a Pohsh space (complete separable 
metric space) [4]: d\{h,k) = Yll^i ^| /o (^(^) ~ ^(■s); ej(s))„(is|, where {ei(s)}^;^ is an 
orthonormal basis for L'^{0,T; Hq). Define 

Am = {heA: h{u:) e Sm, a.s.}. (2.26) 



In order to define the stochastic control equation, we introduce another intensity coef- 
ficient a and also nonlinear feedback forcing R (instead of R) which satisfy 

Condition (C3): (i) a G C([0, T] x V; L{Hq,H)) and there exist constants Kj-i, Ki, 
and Lj, for i = 0,1 and j = 1,2 such that 

\Ht,u)\l^H^^H)<Ko + Ki\u\^+Kn\\u\\Ji, VtG [0,r], Vney, (2.27) 

\a{t,u) -^(i,f)li(j/o,H) < Li\u - v\^ + L2\\u - vf , Vt G [0,r], Vn,uG V, (2.28) 

where \ ■ \l(Ho,h) denotes the (operator) norm in the space L{Hq,H) of all bounded linear 

operators from Hq into H . 

(ii) R : [0,T] X H h^ H is a continuous mapping such that 

\R{t,0)\<Ro, \R{t,u)-R{t,v)\<Ri\u-v\, \/u,v € H, \ft £ [0,T], 

for non-negative constants Rq and Ri. 

Remark 2.3. In contrast with Condition (C2) our hypotheses concerning the control 
intensity coefficient a involve a weaker topology (we deal with the operator norm | • \l(^Ho,h) 
instead of the trace class norm | • \lq)- However we require in (2.27) a stronger bound 
(in the intermediate space 7i). One can see that any noise intensity a satisfies Condition 
(C3)(i) provided Condition (C2) holds with K2 = 0. 

Let M > 0, /i G Am and ^ £ H. Under Conditions (C2) and (C3) we consider the 
nonlinear SPDE with initial condition Uh{0) = ^: 

duhit)+[Auh{t)+B{uhit))+R{t,Uhit))] dt = a{t,Uh{t)) dW{t)+a{t,Uh{t))h{t) dt. (2.29) 

Fix T > and let X := C([0, T]; H^ nL^(0, T; V^ denote the Banach space with the norm 
defined by 

||ti||^ = I sup \u{s)f + \\u{s)fds}\ (2.30) 

Recall that an (.Fi)-predictable stochastic process Uh{t,u)) is called a weak solution in X 
for the stochastic equation (2.29) on [0, T] with initial condition ^ if n S X = C([0, T]; H)r\ 
L^((0, T);y), a.s., and satisfies 

{uh{t),v) - {i,v) + / [{uh{s),Av) + {B{uh{s)) , v) + {R{s,Uh{s)),v)]ds 
Jo 

= {a{s,Uh{s))dW{s),v) + {a{s,Uh{s))h{s) , v) ds, a.s., 
Jo Jo 

for all V G Dom{A) and all t £ [0, T"]. Note that this solution is a strong one in the 
probabilistic meaning, that is written in terms of stochastic integrals with respect to the 
given Brownian motion W. 

The following assertion shows that equation (2.25), as well as (2.29), has a unique 
solution in X, and the X-norm of the solution Uh to (2.29) satisfies a priori bounds 
depending on M when h G Am- 
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Theorem 2.4. Assume that Conditions (C1)-(C2) are satisfied, and that either condi- 
tion (i) or else (ii) below hold true: 

(i) a = a and R satisfies condition (C3)(ii); 

(ii) Condition (C3) is satisfied. 
Then for every M > and T > there exists K2 := K2{T,M) > 0, (which also depends 
on Ki, Ki and Ri, i = 0, 1, and on Kn) such that under conditions IE|^|^ < 00, h £ Am, 
K2 G [0,K2[ and L2 < 2 there exists a weak solution Uh in X of the equation (2.29) 
with initial data Uh{0) = (^ (z H . Furthermore, for this solution there exists a constant 
C := C{Ki, Li, Ki, Kfi, Li, Ri,T,M) such that for h G Am, 

e( sup \uh{t)\^+ [ \\uh{t)fdt+ [ \\uhit)\\j,dt) <C{1 + E\^\^). (2.31) 

^0<t<T Jo Jo ' 

If the constant L2 is small enough, the equation (2.29) has a unique solution in X . If one 
only requires L2 < 2, then equation (2.29) has again a unique solution in X if either a = a 
or if the function h possesses a deterministic bound, i.e., there exists a (deterministic) 
scalar function ip{t) G L^(0, T) such that \h{t)\o < 'ip{t) a.s. 

The proof is similar to that given in [14] for 2D Boussinesq model (2.16) and (2.17) 
(see also [28, 34] for the case of 2D Navier-Stokes equations (2.9)). However, since we deal 
with an abstract hydrodynamical model under a minimal set of hypotheses, the argument 
requires substantial modifications. For the sake of completeness we give a detailed proof 
in an Appendix (see Section 4). Note that only the case a = a will be used for the large 
deviations results. However, the general case is needed in some other frameworks, such as 
the support characterization of the distribution to equation (2.25) in [10]. 

3. Large deviations 

We consider large deviations using a weak convergence approach [3, 4], based on vari- 
ational representations of infinite dimensional Wiener processes. Let e > and let u^ 
denote the solution to the following equation 

du%t) + [Au^t) + B{u%t)) + R{t, u^t))] dt = ^/ecj{t, u^{t)) dW{t) , -^"(0) = ^ G H. (3.1) 

Theorem 2.4 shows that for a any choice of K2 and L2, for e small enough the solution 
of (3.1) exists and is unique in X := C{[0,T],H) n L^{[0,T],V); it is denoted by u" = 
Q^(\feW^ for a Borel measurable function Q"^ : C{[0,T],II) -^ X. A detailed proof of the 
existence of a such function Q^ is given in [32]. 

Let B{X) denote the Borel cj— field of the Polish space X endowed with the metric asso- 
ciated with the norm defined by (2.30). We recall some classical definitions; by convention 
the infimum over an empty set is +00. 

Definition 3.1. The random family {u^) is said to satisfy a large deviation principle on 
X with the good rate function I if the following conditions hold: 

I is a good rate function. The function function I : X ^ [0, cx)] is such that for 
each M E [0, cxo[ the level set {(p £ X : I{(j)) < M} is a compact subset of X. 
For A e B{X), set I {A) = inf„gA/(u). 

Large deviation upper bound. For each closed subset F of X: 

limsup elogP(u^ S F) < -I{F). 

e^O 

Large deviation lower bound. For each open subset G of X: 

lim inf elogP(n= G G) > -I{G). 

£— >0 
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For all h E L'^{[0,T],Ho), let Uh be the solution of the corresponding control equation 
(3.2) with initial condition u/j(0) = ^: 

duh{t) + [Auh{t) + B{uh{t)) + R{t, Uh{t))]dt = a{t, Uh{t))h{t)dt. (3.2) 

Let Co = {/o/i(s)ds : h G L2Qo,T],i/o)} C C([0,r], i^o). Define g^ : C{[0,T],Ho) -^ X 
by G^{g) = "W/i for g = Jq h{s)ds e Co and G^{g) = otherwise. Since the argument below 
requires some information about the difference of the solution at two different times, we 
need an additional assumption about the regularity of the map a{.,u). 

Condition (C4) (Time Holder regularity of a): There exist constants 7 > and C > 
such that for ti, ^2 G [0, T] and u £ V: 

\aituu)-a{t2,u)\LQ <C (l + ||n||)|ti-t2r. 
The following theorem is the main result of this section. 

Theorem 3.2. Suppose the conditions (^Clj and (^C2j with K2 = L2 = are satisfied. 
Suppose furthermore that the conditions (^C3 (ii)J and (C4:) hold. Then the solution (u^) 
to (3.1) satisfies the large deviation principle in X = C{[0,T]; H)n L'^{{0,T)]V), with the 
good rate function 

I,(u)= inf \- f \h(s)\lds]. (3.3) 

{heL^O,T;Hoy.u=gO{f^h{s)ds)}l2 Jo ' ^ '"^ J 

We at first prove the following technical lemma, which studies time increments of the 
solution to a stochastic control problem extending both (3.1) and (3.2). When a, a and 
R satisfy (C2) and (C3), h G Am, the stochastic control problem is defined as in (2.29): 
n|(0) = ^ and 

dul{t) + [Aul{t) + B{ul{t)) + R{t, ul{t))] dt = y/ecr{t, u|(t)) dW{t) + a{t, ul{t)) h{t) dt. 

(3.4) 
To state the lemma mentioned above, we need the following notations. For every integer 
re, let ipn '■ [0,T] — > [0,T] denote a measurable map such that for every s G [0,r], s < 
■0n(s) < (s + c2~") A T for some positive constant c. Given N > 0, h £ Am, and for 
t£ [0,r], let 

GN{t) = [u; : (^supj<(s)(a;)|2) V (^ J \\uUs){u;)fds) < iv}. 

As in Proposition 4.2, we can use a relaxed form of condition (C3 (i)) in order to perform 
calculations in the following lemma; this more general setting is again used in [10]. 

Lemma 3.3. Let eQ,M,N > 0, a satisfy condition (^02^ and a satisfy (2.28) and the 
following condition (3.5) 

\^it,u)\l^Ho,H)<Ko + Ki\u\^ + K2\\u\\\ VtG [0,r], VuGV, (3.5) 

instead of (2.27). Assume that ^ G L'^{Q;H) and let Uh(t) be solution of (3.4) satisfying 
the conclusion of Theorem 2.4- Then there exists a positive constant C (depending on 
Ki,Ki,i = 0,1,2, Lj,Lj,j = 1, 2,Ri,T, M,N,eQ) such that for any h G Am, £ S [Oj^o]; 

/„(/i,e):=E[lG^(r) / |4(s)-<(^„(s))|2dsj <C2-t. (3.6) 

L JO J 

Proof. The proof is close to that of Lemma 4.2 in [14]. However we deal with a class of 
more general functions ipn{s) and do not assume that K2 = L2 = K2 = L2 = 0. As above. 
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to lighten the notation we skip the time dependence of cr, a and R. Let h G Au^ ^ > 0; 
for any s G [0,T], Ito's formula yields 

K(r)-<(s),d4(r))+e / k(4(r))|i^dr. 

J S 

Therefore Inih,e) = I]i<j<6 -^",«' where 



In,l 



2V^E(lG^(r) [ ds f "' {a{uUr))dW{r),uUr)-uUs))), 



In,2 = eE(lG^iT)J dsj W{ul{r))\l^dr 






In,3 = 2e{Ig^^t) J ds J " {a{uUr))h{r),uUr)-uUs))dr), 
/„,4 = -2E(1g^^t)J dsj (AuUr) , ulir) - ulis)) drj, 
In,, = -2e{Ig,^t) I ds r\B{ul{T)),ul{r)-ul{s))dr), 

Infi = -2E(lG^^T)jdsj^^\R{ul{r)),ul{r)-ul{s))dr). 

Clearly Gn{T) C GN{r) for r £ [0,T]. In particular this means that |u^(r)| + |u^(s)| < N 
on GN{r) for < s < r < T. We use this observation in the considerations below. 
The Burkholder-Davis-Gundy inequality and (C2) yield for < e < eo 

l^,i| < 6^^^ d.E(^ "' |a«(r))|i^lG^(,)|<(r)-4(s)|2dr)" 

< 6y^2^ f dsE(^f"\Ko + K,\ul{r)\^ + K2\\ul{r)f]dry. 

Schwarz's inequality and Fubini's theorem as well as (2.31), which holds uniformly in 
e G]0,eo] for fixed Eq > since the constants Ki and Li are multiplied by at most eq, 
imply 



|/n,i| < 6^2eoiVr [e / {Ko + K^ \ul{r)\^ + K2 \\ul{r)f) ( [ ds) dr 

L Jo ^ J{r-c2-")y0 ' 

< Ci2-^ (3.7) 

for some constant Ci depending only on Ki, Ki,i = 0, 1, 2, Lj, Lj,j = 1, 2, Ri, M, eq, N 
and T. The property (C2) and Fubini's theorem imply that for < e < eo, 

\In,2\ < EE(lG^f^T)jdsj^\Ko + K^\ul{r)\'' + K2\\ul{T)f)dr) 

< EoE [ lG^^T){Ko + KiN + K2\\uUr)f)c2-''dr < C22~'' (3.8) 
Jo 

for some constant C2 depending on Ki, i = 0,1,2, eq, N and T. Schwarz's inequality, 
Fubini's theorem, (C2) and the definition (2.26) of Am yield 



|I„,3|<2E 1g^(T) / ds 



T 
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XI {Ko + K,\ui{r)f + K2\\uUr)fy' \hir)\o\ ui{r) - <(s)| dr 



i: 

Jo 



lGr,iT)\h{rMKo+KiN + K2\\ul{r)r)-2{ I ds] dr 

^ J(r-c2-")V0 



< 4ViVc2^"VME 1 



LG^(T) I (Ko + KiN + K^WuUrW) dr) < G, 2"", (3.9) 

for some constant C3 depending on Ki, i = 0, 1,2, M and A^. Using Schwarz's inequality 
we deduce that 

rT ci'nM 



In A < 2E(lG^(r) / ds 



Js 

T 



< ^e(1g^^t)I ds\\ul{s)f j 



dr[-\\ul{r)r +\\ul{r)\\\\ul{s)\\\\\ 
dr] <ciV2-("+^). 



(3.10) 



The antisymmetry relation (2.2) and inequality (2.7) yields 

\{B{uUr)),ulir)-uUs))\ = |(S«(r)) , 4,(.)>| < ||<(r)f + C7|<(r)|2||<(.)||^. 
Therefore, 



In.5\ <2E(1 



LGjv(T) 



ds / dr ||<(r)f 

is 

T 



r-tZ-nCs) 



ef^\l2\ ^ r(l) , r(2) 
- -'n,5 "I" -"n.S- 



ds 



'{r-c2-")V0 



+ 2C¥.(lG,(T)j ds\\ul{s)tnj^ " V|<(r)| 
Fubini's theorem implies 

/^ < 2E(la^^r)rdr\\uiir)f f 
^ Jo J{i 

< 2c2-^¥.(lG,(T)j dr ||<(r)f ) < CN2-^. 

Using (2.3), we deduce that on Gn{T) we have 

\\uUs)\\tids < al sup \uUs)f f\\nUs)fds<alN'. 
se[o,T] Jo 



(3.11) 



(3.12) 



Thus 



li^l < 2CNe(^Ig^^t)J^ ds||4(s)||^)c2-" < 2alCN^c2-^. (3.13) 



Finally, Schwarz's inequality implies that 



Infil <2E 



1g^{T) / ^s^ "\Ro + Ri\ul{r)\){\ul{r)\ + \ul{s)\)dr 



<CN^2 



2r,—n 



(3.14) 



D 



Collecting the upper estimates from (3.7)-(3.14), we conclude the proof of (3.6). 

In the setting of large deviations, we will use Lemma 3.3 in the case when then a = a 
satisfies Condition (C2) with K2 = L2 = and with the following choice of the function 
ipn- For any integer n define a step function s 1— > s^ on [0, T] by the formula 

Sn = tk+i = {k + 1)T2-'' for sE [A:r2-",(A; + l)r2~"[. (3.15) 

Then the map ipn{s) = Sn clearly satisfies the previous requirements with c = T. 
Now we return to the setting of Theorem 3.2. 
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Let £0 > 0, (/le, < e < Eq) be a family of random elements taking values in the set Am 
given by (2.26). Let u^^, or strictly speaking, u| , be the solution of the corresponding 
stochastic control equation with initial condition Uh^ (0) = ^ E H: 

duh^ + [Auh^ + B{uh,) + R{t,uyJ]dt = a{t,Uh^)h,{t)dt + ^ a{t,Uh,)dW{t). (3.16) 

Note that for W^ = W, + -^ J^^he{s)ds we have that Uh^ = Q''{^/eW^). The following 
proposition establishes the weak convergence of the family {uh^) as e — > 0. Its proof is 
similar to that of Proposition 4.3 in [14], but allows time dependent coefficients R and a. 

Proposition 3.4. Suppose that the conditions (CS-) and (^C2j are satisfied with K2 = 
L2 = 0. Suppose furthermore that R and a satisfy the conditions (C3)(ii) and (C4). Let 
^ he J^Q-measurahle such that E\S\'^^ < +00, and let /i^ converge to h in distribution as 
random elements taking values in Am, where this set is defined by (2.26) and endowed 
with the weak topology of the space ^2(0, T; Hq). Then as e ^ 0, the solution Uh^ of (3.16) 
converges in distribution to the solution u^ of (3.2) in X = C{[0,T]; H) H L'^{{0,T);V) 

endowed with the norm (2.30). That is, as e ^ 0, ^^( ^/£{W, -\- -^ J^ h^{s)ds^ j converges 

in distribution to ^°( J^ h{s)ds) in X. 

Proof. Since Am is a Polish space (complete separable metric space), by the Skorokhod 
representation theorem, we can construct processes {h£,h,W^) such that the joint dis- 
tribution of {h£,W^) is the same as that of {hi;,W^), the distribution of h coincides 
with that of h, and /i^ -^ h, a.s., in the (weak) topology of Sm- Hence a.s. for every 
t £ [0,T], Jq hi;{s)ds — Jq h{s)ds — > weakly in Hq. To lighten notations, we will write 
{h„h,W') = {h„h,W). 

Let Ue = Uh^ — Uh] then Us{0) = and 

dUe + [AUe + B{uh^) - B{uh) + R{t,Uh^) - R{t,Uh)]dt 

= [a{t,Uh,)he - a{t,Uh)h]dt + ^/e a{t,Uh,)dW{t). (3.17) 

On any finite time interval [0, t] with t <T, Ito's formula, (2.8) with ?? = 2 ^^d condition 
(C2) yield for e > 0: 

\Us{t)\^ + 2 [ \\Ue{s)fds = -2 [ {B{uhM) - B{uh{s)) , Ue{s))ds 
Jo Jo 

- 2 / {R{S, Uh, (s)) - Ris, Uhis)) , Ue{s)) ds 

Jo 

+ 2 / {a{s,UhSs))h£{s) - a{s,Uh{s))h{s),Ue{s))ds 
Jo 

+ 2^8 [ {Us{s),a{s,Uh,is))dW{s))+e [ \ais,UhM)\lnds 
Jo Jo 

r-t 3 »i 

/ \\Ue{s)f ds + y^Ti{t,e) + 2 / [Ci\\uh{s)tn + Ri + ^/ri\h,{s)\o)\Ue{s)\^ds, 
Jo .^1 Jo 2 



< 



where 



ri(t,e) = 2^8 f {Ue{s),a[s,UhAs))dW{s)), 
Jo 

T2{t,e) = e [ iKo + Ki\uhM\^)ds, 
Jo 
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ft 



r3(t, e) = 2 y (ais, Uh{s)) {Kis) - his)) , Ueis)) ds. 
This yields the following inequality 

\Ueit)f+ / \\Ue{s)fds<y^Ti{t,e) + 2 / Ci\\uh{s)\\^+Ri + y^i\he{s)\o \Ue{s)\Us. 
Jo I^ -^0 L 2 J 

(3.18) 
We want to show that as e — > 0, ||f4||x — > in probability, which implies that Uh^ — > Uh 
in distribution in X. Fix A^ > and for t £ [0,T] let 

GNit) = I sup \uh{s)\^ < n} n \ [ \\uh{s)fds<N}, 

^ 0<s<t ^ ^ Jo ' 

GN,e{t) = GN{t)n{ snp \uhM\^<N}n\ [ \\uhMfds<N}. 

^ 0<s<t ^ ^ Jo ^ 

The proof consists in two steps. 
Step 1: For any eq g]0, 1]0, sup sup F^GnAT)! ^ as N ^ oo. 

Indeed, for e £]0,eo], h,hs S Am, the Markov inequality and the a priori estimate (2.31), 
which holds uniformly in e G]0,eo], imply 

F{GN,e{TT) < ^( sup \uh{s)\^ > A^) +P( sup \uhM^ > n) 

^0<s<T ' ^0<s<T ' 

+ P(y" \\uh[s)fds>N)^F^ ||nft^(s)||2)ds>iv) 
<— sup e( sup \uh{s)\^ + sup \uhAs)\^+ {\\Ms)f + \WhAs)f)ds) 

^^ h,he&AM ^0<s<T 0<s<T Jo ^ 

<C(l + E|^|^)7V-\ (3.19) 

for some constant C depending on T and M. 

Step 2: Fix N > 0, h,h^ £ Am such that as e ^ 0, /i^ ^ /i a.s. in the weak topology of 

L^(0, T; Hq); then one has as e ^ 0: 

rT 



E 



0. 



(3.20) 



1g^,(T) sup \UM'+ / WeitWdt 

^ 0<t<T Jo 

Indeed, (3.18) and Gronwall's lemma imply that on GN,e{T), 
sup \Ue{t)\'^ < \ sup (ri(t,e)+r3(t,e)) +eCJ exp (2aoCi A^^ + 2RiT + 2^J LiMT) , 



0<t<T 



-0<t<T 



where G^ = T{Kq + KiN). We also use here the fact that by (2.3) 

T r-T 

\\uh{s)\\'^ds < aQ sup \uh{s)\^ / \\uh{s)\\^ ds < aoN"^ oiiGN,e{T). 
se[o,T] Jo 

Using again (3.18) we deduce that for some constant G = G(T, M, N), one has for every 
£>0: 



E(l 



GnAT) 



\Usfx) <c(e + E 



GnAT) 



sup {n{t,e)+n{t,e)) 



0<t<T 



(3.21) 



Since the sets Gn^£{.) decrease, E(1g'^^('p) supQ<j<y |Ti(t,e)|) < E(Ae), where 
Ae:=2Vesup / lGj^^(s)(Ue{s),a{s,Uh,{s))dW{s) 

0<t<T Jo ' ^ 
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The scalar-valued random variables Ag converge to in L^ as e — > 0. Indeed, by the 
Burkholder-Davis-Gundy inequality, (C2) and the definition of Gn^s{s), we have 



E(A,) < 6ViE{ / 1g^^^(,) \U,{s)\^ |a(s,n;,,(s))|i^ds} 



<6V^e[[4nI lG^^^is){Ko + Ki\uhAs)f)dsy]<C{T,N)^. (3.22) 

In further estimates we use Lemma 3.3 with ?/^„ = s„, where Sn is the step function 
defined in (3.15). For any n,N > 1, if we set tk = kT2~^ for < A; < 2", we obviously 
have 



0<t<T 



i=l 



e(1g^,(t) sup |r3(t,e)|) <2 Vf,(Ar,n,e) + 2E(t5(iV,n,e)), (3.23) 

where 
ri(iV,n,e) =E 



lG^,.(r) sup 
o<t<r 



O" 



{s,Uh{s)){K{s) - his)) , [[/,(s) - Ue{Sn)])ds 



T2(iV,n,e)=E 



T3(iV,n,e)=E 



r4(7V,n,e) =E 



GN,e{T) 



X sup 

0<t<T 



[a{s,Uhis)) -a{sn,Uh{s))]{he{s) -h{s)), Ueisn))ds 



Gn,AT) 



X sup 

0<t<T 



[a{Sn,Uh{s)) - a{Sn,Uh{Sn))] {he{s) - h{s)) , C4(s„) jds 



'^GN^iT) sup sup 

l<fc<2" tfc_i<t<tfe 



0-(ifc, Wft(tfc)) / (/ie(s) - h{s)) ds , [/^(tfc) 



r5(7V,n,e) =lG^^^(T)Y,\((^{tk,Mtk)) / (/ie(s) - /i(s)) ds,Ueitk) 



k=l 



Using Schwarz's inequality, (C2) and Lemma 3.3 with ipn = Sn, we deduce that for some 
constant Ci := C{T,M, N) and any e E]0,eo], 

fi{N,n,e) < Ehc^^^iT) I {Ko + Ki\uhis)\^)^^ \he{s) - h{s)\o \Ue{s) - Ue{sn)\ds 



< E 



IGat.CT) / {\uhAs) - '^h,{Sn)\'^ + \uhis) - Uh{Sn)\'^} ds 

Jo 



T 



X \/2(Ko + KiN) f E / |/i£(s)-/i(s)|gdsj' <C7i 2-t. 



1 



(3.24) 



A similar computation based on (C2) and Lemma 3.3 yields for some constant C3 := 
C(r, M, A^) and any e g]0, eo] 



f3(iV,n,e)< V2A^Li(e iG^^(r) / |n;,(s) - n;,(s„)|2 ds ) Ve / |/i,(s) - /i(s)|2ds 



1 , /-T 

2 



|2j„\2 



<C32-?. (3.25) 

The Holder regularity (C4) imposed on cr{.,u) and Schwarz's inequality imply that 



T 



T2{N,n,e) <CVN2~'^^E[lG^^^iT) J {1 + \\uhis)\\)\he{s)-h{s)\dsj <C22-'^^ (3.26) 
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for some constant C2 = C{T, M,N). Using Schwarz's inequality and (C2) we deduce for 

C4 = C{N, M) and any e g]0, Eq] 



fi{N,n,e) <E 



Ig^JT) sup {KQ + Ki\uh{tk)?)'^ r \he{s)-h{s)\Qds\Ue{tk)\ 

l<fe<2" Jtk-1 

'■tk 



<2ViV(A'o + i^iiV)]E( sup / |/ie(s)-/i(s)|o(is) <4(74 2-t. (3.27) 

Finally, note that the weak convergence of h^ to h implies that for any a, 5 G [0, T], a < 6, 
the integral / hi;(s)ds -^ J h{s)ds in the weak topology of Hq. Therefore, since for the 
operator a{tk,Uh{tk)) is compact from Hq to H, we deduce that for every k, 



a{tk,Uhitk)) 



hs{s)ds 



h{s)ds 



as e 



0. 



H 



Hence a.s. for fixed n as e — > 0, T^{N,n,e,uj) -^ 0. Furthermore, T^{N,n,e,uj) < 
C{Kq, Ki, N, M) and hence the dominated convergence theorem proves that for any fixed 
n,N, ¥.{f^{N,n,£)) ^ as e ^ 0. 

Thus, (3.23)-(3.27) imply that for any fixed A^ > 1 and any integer n > 1 



lim sup E 



1gjv,(t) sup \T'i{t,e)\ 

0<t<T 



<Cntm2-<^''\\ 



Since n is arbitrary, this yields for any integer A^ > 1: 



limE 

e-»0 



^GnAT) sup |T3(t,e) 
o<t<r 



0. 



Therefore from (3.21) and (3.22) we obtain (3.20). By the Markov inequality 



\U, 



for any (5 > 0. 



Finally, (3.19) and (3.20) yield that for any integer A^ > 1, 

limsupP(||[/e||x >6) < C{T,M)N-^, 

for some constant C{T,M) which does not depend on N. This implies linig^oIPdlt^llx > 
(5) = for any 6 > 0, which concludes the proof of the proposition. D 

The following compactness result is the second ingredient which allows to transfer the 
LDP from \feW to u^ . Its proof is similar to that of Proposition 3.4 and easier; it will be 
sketched (see also [14], Proposition 4.4). 

Proposition 3.5. Suppose that (^ClJ and (^C2j hold with K2 = L2 = and that conditions 
(C3)(ii) and (C4) hold. Fix M > 0, ( € H and let Km = {uh £ X : h £ Sm}, where u^ 
is the unique solution of the deterministic control equation (3.2) and X = C{[0,T];H) f] 
L^(0, T;y). Then Km is a compact subset of X. 

Proof. Let {un} be a sequence in Km, corresponding to solutions of (3.2) with controls 
{hn} in Sm- 

dUn{t) + [AUn{t) + B{Un{t)) + R{t, Unit))] dt = a{t, Un{t))hn{t)dt, U„(0) = ^. 

Since Sm is a bounded closed subset in the Hilbert space L^(0, T; Hq), it is weakly compact. 
So there exists a subsequence oi {hn}, still denoted as {hn}, which converges weakly to a 
limit h in -L^(0, T; Hq). Note that in fact h G Sm as Sm is closed. We now show that the 
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corresponding subsequence of solutions, still denoted as {u„}, converges in X to li which 
is the solution of the following "limit" equation 

du{t) + [Au{t) + B{u{t)) + R{t, u{t))]dt = a{t, u{t))h{t)dt, u{0) = ^. 

This will complete the proof of the compactness of Km- To ease notation we will often 
drop the time parameters s, t, ... in the equations and integrals. 

Let Un = Un — u; using (2.8) with t] = ^, (C2) and Young's inequality, we deduce that 
forte [o,r], 

\Un{t)\^ + 2 f\\Un{s)fds = 

Jo 

rt i-t 

2 / {B{Un{s))-B{u{s)),Un{s))ds-2 / {R{s,Un{s)) - R{s,Uh{s)),Un{s))ds 



+ 2 



I y [cr(s, Un{s)) - Cr(s, u{s))\ /l„(s), C/„(s) 

+ (CJ(S, U{s)){hn{s) - h{s)) , f/„(s)) }ds 
< [ \\Un{s)fds + 2 I \Un{s)WC.\\u{s)\\\^+Ri + ^i\K{s%)ds 

Jo Jo ^ 

+ 2 / (cj(s, u{s)) [hn{s) - h{s)] , t/„(s)) ds. (3.28) 



The inequality (2.31) implies that there exists a finite positive constant C such that 



sup 

n 



sup (|n(t)|2 + |u„(t)|2)+ / {\\u{s)f + \\u{s)f^ + \\ur,{s)f)ds =C. (3.29) 

0<t<T JO -J 



Thus Gronwall's lemma implies that 

/•T 5 

snp\Unit)\^+ / \\Un{t)fdt<exp(2{CiC + RiT+y^LjIf)) Y^iIm. 

t<T Jo \ 2 / ^— ' 



(3.30) 



i=l 



where, as in the proof of Proposition 3.4, we have: 



L 



1 

n,N 



o- 







{s,u{s)) [hn{s) - h{s)] , Un{s) - UnisN))\ds, 



^n,N 



'n,iV 



'n,iV 



T 



T 



[a{s, u{s)) - cr{sN, u{s))] [K{s) - h{s)] , Un{sN 

[a{sN,u{s)) - (7{sn,u{sn))] [hn{s) - h{s)] , Un{sN 



ds, 



ds, 



sup sup 

l<fc<2^tfe_i<i<ifc 



o-{tk,u{tk)) I {he{s) - h{s))ds , Un{tk] 






•-n,N - '^{(^iik,u{tk)) [hn{s)-h{s)]ds ,Un{tk)y 

k=l -^^fc-i 

Schwarz's inequality, (C2) and Lemma 3.3 imply that for some constants Cj, which do 
not depend on n and A^, 

/^,Ar<Co( / \hn{s) - h{s)\ldsy i / [\Un{s) -Un{sN)\'^ + \u{s)-u{sN)\'^)ds 



^JO 

N 

< Ci 2~- , 



(3.31) 
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rT . 1 , /-T , 1 

\u{s) -u{sN)\'^dsy y / 

In.N < C'o 



In,N<Co[l \u{s)-u{sN)\^dsyU |/i„(s) - /i(s)|g ds) ' < Cg 2^f , (3.32) 

l+( sup \u{t)\) sup (|u(t)| + |n„(t)|)l 2~f < C4 2~f . (3.33) 

Furthermore, condition (C4) implies that 



/nW<Co2-^T sup {\u{t)\ + \Un{t)\) (1 + ||n(s) ||)(|/l(s)|o + |/ln(s)|o) ds < ^2 2"^^. 

0<t<T Jo 

(3.34) 
For fixed N and k = 1,- ■ ■ , 2^, as n ^ oo, the weak convergence of /i„ to /i implies that of 
/j* {hn{s) — h{s))ds to weakly in Hq. Since a{u{ti^)) is a compact operator, we deduce 

that for fixed k the sequence a{u{tk)) f^ ^ {hn{s) — h{s))ds converges to strongly in H 

as n — > oo. Since sup„ ^ \Un{tk)\ < 2vC, we have lim„/^^ = 0. Thus (3.30)-(3.34) yield 
for every integer iV > 1 

limsup \ sup \Un{t)\^ + I \\Un{t)f dt \ < C2"^(i'^3). 

n->oo \^t<T Jo J 

Since N is arbitrary, we deduce that ||C/n|LY -^ as n — > oo. This shows that every 
sequence in Km has a convergent subsequence. Hence Km is a sequentially relatively 
compact subset of X. Finally, let {un} be a sequence of elements of Km which converges 
to V in X. The above argument shows that there exists a subsequence {un^, k > 1} which 
converges to some element u^ G Km for the same topology of X. Hence v = u^, Km is a 
closed subset of X, and this completes the proof of the proposition. D 

Proof of Theorem 3.2: Propositions 3.5 and 3.4 imply that the family {u^} satisfies 
the Laplace principle, which is equivalent to the large deviation principle, in X with the 
rate function defined by (3.3); see Theorem 4.4 in [3] or Theorem 5 in [4]. This concludes 
the proof of Theorem 3.2. □ 

4. Appendix : Proof of the Well posedness and apriori bounds 

The aim of this section is to prove Theorem 2.4. We at first suppose that conditions 
(C1)-(C3) are satisfied. Let F : [0,T] x V ^ V be defined by 

F{t, u) = -Au - B{u, u) - R{t, u) , yte [0, T], Vn G V. 

To lighten notations, we suppress the dependence of a, a, R and F on t. The inequality 
(2.8) implies that any 77 > there exists C,, > such that for u,v ^V , 

{F{u) - F{v) ,u-v)<-{l- r])\\u - vf + {Ri + CJv\\j^) \u - v\'^. (4.1) 

Let {(pn}n>i be an orthonormal basis of the Hilbert space H such that (pn G Dom{A). 
For any n > 1, let Hn = span{ipi,- ■ ■ ,(pn) C Dom{A) and Pn '■ H ^> Hn denote the 
orthogonal projection from H onto Hn, and finally let an = Pn<y and ct„ = Pn<5"- Since P„ 
is a contraction oi H, from (2.22) we deduce that \cyn{u)\\„ < l'^('")lio- 

For h £ Am^ consider the following stochastic ordinary differential equation on the 
n-dimensional space Hn defined by ttn,h(0) = Pn^-, and for v G Hn- 

d{Un,h{t),v) = [{F{Un,h{t)),v) + {d{Un,h{t))h{t), v)] dt + {a{Un,h{t))dWn{t),v), (4.2) 

where Wn{t) is defined in (2.23). Then for A; = 1, • • • , n we have 

d{Un,h(t),ipk) = [{F{Un,h{t)), ^k) + {a{Un,h{t))h(t), ipk)] dt 
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+ ^Qj {(^iun,h(t))ej , ^k)dl3j{t). 

Note that for v G Hn the map u G Hn ^^ {Au + R{u) , v) is globally Lipschitz uniformly 
in t, while using (2.7) we deduce that the map u G Hn i— > {B{u) , v) is locally Lipschitz. 
Furthermore, since there exists some constant C{n) such that \\v\\ < C{n)\v\ for v G Hn, 
Conditions (CI) and (C2) imply that the map u G Hn i— > [{an{u)ej , (pk) '■ I < j,k < n), 
respectively u G Hn ^^ [{an{u)h{t) , tpk) : 1 < k < nj, is globally Lipschitz from Hn to 
n X n matrices, respectively to M" uniformly in t. Hence by a well-known result about 
existence and uniqueness of solutions to stochastic differential equations [22], there exists 
a maximal solution Un,h = Ylk=ii'^n,h , ^k) 'fk to (4.2), i.e., a stopping time Tn^h < T such 
that (4.2) holds for t < Tn,h and as t t ^11,11 < T, \un,h{t)\ — > 00. 

The following proposition shows that Tn,h = T a.s. It gives estimates on Un,h depending 
only on T, M, Ki, Li and E|,^p'P, which are valid for all n and all K2 G [0,-^2] for some 
K2 > 0. Its proof depends on the following version of Gronwall's lemma (see [14], Lemma 
3.9 for the proof). 

Lemma 4.1. Let X, Y , I and (p he non-negative processes and Z he a non-negative 
integrahle random variable. Assume that I is non- decreasing and there exist non-negative 
constants C , a, (5, 7, 5 with the following properties 

i-T 

/ Lp{s) ds<C a.s., 2l3e^ < 1, 25e'^ < a, (4.3) 

Jo 

and such that for < t <T, 

X{t) + aY{t) < Z+ [ <f{r)X{r)dr + I{t), a.s., 

Jo 

E(/(t)) < PK{X{t))+j [ K{X{s))ds + 6E{Y{t))+C, 

Jo 

where C > is a constant. If X & L°°([0,T] x il.), then we have 

K[X{t) + aY{t)] <2e^p{C + 2tje^){E{Z) + C), tG[0,r]. (4.4) 

The following proposition provides the (global) existence and uniqueness of approximate 
solutions and also their uniform (a priori) estimates. This is the main preliminary step in 
the proof of Theorem 2.4. As in Lemma 3.3 we can made this step under less restrictive 
growth conditions concerning a than (2.27) in (C3). 

Proposition 4.2. Fix M > 0, T > and let Conditions (^Clj-(^C3j he in force with the 
assumption (3.5) instead of (2.27). For any integer p > 1 there exists K2 = K2{p,T,M) 
(which also depends on Ki, Ki and Ri, i = 0, 1, and on K2 if K2 7^ K2) such that the 
following result holds. Let h G Am, < K2 < K2 and ^ G L'^^{Q.,H). Then equation (4.2) 
has a unique solution on [0, T] (i.e. Tn,h = T a.s.) with a modification Un,h £ C{[0,T],Hn) 
and satisfying 

supe( sup lunMt)]^" + [ \\nn,h{^)fWnAs)f^'"'^ds)<C{mfP + l) (4.5) 

n V 0<t<T Jo ' 

for some positive constant C (depending on p, Ki,Ki,i = 0, 1, 2, Rj,j = 0, 1, T, M). 

Proof. Let Un^t) be the approximate maximal solution to (4.2) described above. For 
every A^ > 0, set 

T^=inf{i: |'U„,,,(t)| >iV}Ar. 
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Let Hn : Hq ^ Hq denote the projection operator defined by UnU = X]fc=i (^' ^fc) Cfc, 
where {e^, fc > 1} is the orthonormal basis of H made by eigen-elements of the covariance 
operator Q and used in (2.23). 

Ito's formula and the antisymmetry relation in (2.2) yield that for t G [0, T], 

/■tArjv /"tATjv 

|u„,h(tATJv)|^ = |Pn^|^ + 2 / {aniUn,his))dWnis),Un,his)) - 2 / ||n„,/,(s) f ds 

Jo Jo 

-2 / {RiUn,his)) - ^n{Un,his))his),Un,hi^)) ds + / |<7n(lir!,,/i(s)) n^li^ ds. 

(4.6) 

Apply again Ito's formula for f{z) = z^ when p > 2 and z = \un,h{'t A T7v)p. This yields 
for t G [0,T], and any integer p > 1 (using the convention p{p — l)x^^^ = if p = 1) 

|n„,;.(tAr;v)P^ + 2p/ h„./.(r)|2(P-i) ||n„.,,(r)f dr < |P„C|2P + V r,(t), (4.7) 
-'0 ,77:.^ 



l<j<5 



where 



Ti{t) = 2p / {Ro + Ri\un,h{r)\)\un,h{r)\^P-^dr, 

Jo 

/■tATjv 

T2it) = 2p / {an{unMr))dWn{r),Un,h{r))\un,h{r)\^^''-'^ 

Jo 

/■iATjv 

r3(t) = 2p |(^n(n„,„(r))Mr),n„.;,(r))||n„,^(r)|2(P-i)dr, 

JO 



r4(t) 



p 



tArjv 



|a„(n„,;,(r)) n„|2 |u„_;,(r)|2(P-i)dr, 



r5(t) = 2p(p-l) / \nr,aUnn,h{r))urr,h{r)\lWn,hir)\^^''^^dr. 

Jo 

Since h S Am, the Cauchy-Schwarz inequality and condition (3.5) imply that 

T^{t) < 2p I NKo + JKi K,h{r)\ + yjK2 \\un,h{^)\\] |/i(r)|o \un,h{r)\^"'^dr 



< P 





tArjv 



\un,h{r)fK.hir)\^^'"^^dr + 2pK2 / \h{r)\l\un,h{r)\^P dr 

Jo 



2 JO 

I rtATff I f'tf^TN 

+ 2pJki / \h{r)\^\un.h{r)\''''dr + 2pJkQ / \h{r)\Q\un,h{r)?''~^dr. 
Jo Jo 

Therefore using the inequality |npP~^ <\ + \u\'^^to bound the last term we obtain 



T3{t) < 



p 



\un,h{r)f\unMr)\^'-''-^Ur + 2pJko / \h{r)\odr 



+ 2p 



tArjv 



ko + Jki]\h{r)\o+k2\h{r] 



\u. 



,Mr)\'''dr. (4.8) 



Using condition (C2), relation (2.22) and also the fact that 

hiu)\\c{Ho,H) = h*{u)\\c{H,Ho) < Wiu)\LQ, 

we deduce that 

/■tArjv 

T4it)+nit) < {2p^-p)K2 / \\un,h{r)f\un,h{r)?^-^^dr 

Jo 
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+ {2p'-p) J [K^\u.,,^,{r)\'P + Ko\un,h{r)\"^''-'^)dr 

rtATN 

<{2p^-p)K2 / \\nn.h{r)f\un,h{r)\"^'"'^dr 

Jo 



+ Cb 



t/\TlsI 



[ifo + (i^o + i^i)hn,h,(r)|2P] dr. 



Consequently (4.7) for K2 < (4p - 2)"^ yields 

/■tATjv 

Jo 



(4.9) 



{Ko + Ro)T+JKo / |/i(r)|odr 



< l^n^P^ + Cp 



for i G [0,T], where I{t) = supo<s<j 1^2(8)1 and 



tArjv 



+ / ip{r)\un,hir)\^Pdr + I{t) 
'0 



^{r) =Cp{Ro + Ri + Ko + Ki + 



'Ko + \^Ki 



\hir)\o + K2\h{r)\l 



for some constant Cp > 0. The Burkholder-Davies-Gundy inequality, (C2) and Schwarz's 
inequality yield that for t £ [0, T] and /? > 0, 



EI(t) =e( s^upjr2(s)|) < 6pE[f ^"^ \un,h{r)\^^^'"^^ Wn{un,h{r)) Un\l^ drj 
<Pk( sup |n„,,(5)P^) + ^^ E r^" |k„,/^(r)f |n„,,(r)|2(P-i)dr 

^0<s<tATjv ^ P JO 



V(Ko+i^l 



-E 



tArjv 



Thus we can apply Lemma 4.1 for 



\un,hir)\^^dr + 



9p^Ko 

/9 



T. 



(4.10) 



X(t)= sup \Un,h{s)\^^, Y{t)-- 

0<s<tATN JO 



tATjv 



„,(r)||2|n„,,(r)p(P-i)dr. (4.11) 



All inequalities for the parameters (see (4.3)) can be achieved by choosing K2 small enough. 
Thus there exists K2 such that for < K2 < K2 we have 

supE( sup K^hl^P + r \\un,h{s)\\ |n„,/.(s)p(P-i) ds) < C{p) 

n ^0<s<rjv Jo ^ 

for all n and p, where the constant C(p) is independent of n. 

Now we are in position to conclude the proof of Proposition 4.2. As A^ ^ cxo, tat | t„^/j, 
and on the set {t^.h < T}, we have supo<^<^^ \un,his)\ -^ 00. Hence P(t„ j,, < T) = and 
for almost all u, for A^(a;) large enough, T]^(^^^{u)) = T and m„^/j(.)(u;) G C([0, T], i/„). By 
the Lebesgue monotone convergence theorem, we complete the proof. D 

Remark 4.3. If K2 < 2 — 6 for some 5 > then the bound K2 does not depend on 
K2. Indeed, one may slightly change the proof by replacing the inequality in (4.8) by the 
following 

rtATN 

Ut) <pK2 \\un,h{r)f kn,h(r)p(^-^) dr 

Jo 

i-tArpf 



+ C„ 







'Ko\h{r)\o + 



A'o + V^i 



\h{r)\o + \h{r)\l] \u.,,,{r)\'P 



dr. 
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Therefore we can exclude the first term of right hand side and obtain relation (4.9) with 
coefficients independent of K2 ■ 

Proof of Theorem 2.4: 

Let us at first suppose that condition (C3) is satisfied. Let $7^ = [0, T] x $7 be endowed 
with the product measure ds (8) dF on ;B([0, T]) J^. Let K2 be defined by Proposition 
4.2 with p = 2. The inequalities (4.5) and (2.3) imply that for K2 G [0, K2] we have the 
following additional a priori estimate 

rT 



supE / \\u,,^his)\\^ds < C2{1 + E|e|^). (4.12) 

n Jo 



'0 

The proof consists of several steps. 

Step 1: The inequalities (4.5) and (4.12) imply the existence of a subsequence of 

iun,h)n>o (still denoted by the same notation), of processes 

UheX := L^i^T, V) n L^{QT,n) n L^{Q, i°°([0, T],H)), 

Fh G L'^{yiT,V') and Sh,Sh G L'^{Qt,Lq), and finally of random variables Uh{T) G 

L^(0,i?), for which the following properties hold: 

(i) Un,h -^ Uh weakly in L'^{Qt, V), 

(ii) Un,h -^ Uh weakly in L^{Qt,'H), 

(iii) Un,h is weak star converging to u^ in L^{Q, L°°{[0,T], H)), 

(iv) Un,h{T) -^ Uh{T) weakly in L'^{Q,H), 

(v) F{un,h) -^ Fh weakly in L^i^r, V), 

(vi) cFn{un,h)^n ^ Sh Weakly in L'^{Q.t,Lq), 

(vii) dn{un,h)h — > Sh weakly in L^{ViT-,H). 

Indeed, (i)-(iv) are straightforward consequences of Proposition 4.2, of (4.12), and of 
uniqueness of the limit of E J„ (un^h{'t),v{t))dt for appropriate v. Furthermore, given v G 
L'^{nT,V), we have Av G L'^{nT,V'). Since for u,v e L^{nT,V), E J^ {Au{t) , v{t)) dt = 
E/(f(u(i), Av{t))dt, 

E [ {Aun,h(t),v{t))dt^ E [ {Auhit) , vit)) dt. (4.13) 

Jo Jo 

Using (4.5) withp = 2, (2.4), (4.12), condition (C3), the Poincare and the Cauchy-Schwarz 
inequalities, we deduce 

SUpE / \{B{Un,h{t)),<i)) + iRiUn,hit)),v{t))\dt 
n Jo 

T 

2\ 



<CiSUp|e/ \\Un,hmti+E [ \Un,h{t)\^dt]+ C2E [ {l + \\v{t)\\'')dt 
n I Jo Jo ) Jo 

< Cafl + ^lel^ + E / ||t;(t)fdt^ <+oo. 



Hence {B{un,h{t))+R{un,h{'t)) , n> 1} has a subsequence converging weakly in L^{Q.t, V'), 
which completes the proof of (v) . 

Since n„ contracts the | • | norm, (C2), (2.22) and (4.5) for p = 2 prove that (vi) is a 
straightforward of the following 

supE / k„(w„,/,(t))n„|i dt < iror + supE / (Ki|n„,/,(t)p + K2\\un^h{t)f) dt < 00. 

n Jo n Jo 
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Finally, using (2.27) in (C3), Holder's inequality, (4.5) with p = 2 and (4.12), we deduce 

f^ 4 fT I— I— rZ^ 4 4 

E / |<5-„(M„,/,(s)/i(s)|3(is <E / \\jKQ^\]Kx\u^.h{,s)\^^Kn\un^h\-H\^ \h[s)\l ds 
Jo Jo 

<Ci(eJ \h{s)\ldsy(EJ [l + \un,hist + \K,hmn]dsy <C{M,T) 

for every integer n > 1. This completes the proof of (vii) . 

Step 2: For 5 > 0, let f £ H^{-5,T + 6) be such that ||/||oo = 1, /(O) = 1 and for 
any integer j > 1 set gj{t) = f{t)(pj, where {^j}j>i is the previously chosen orthonormal 
basis for H. The Ito formula implies that for any j > 1, and for < t < T, 

4 
i=l 

where 

rT rT 

^1,3= I iUn,h{s),fj)f'{s)ds, ll,j= I {(yn{Un,h{s))TlndW{s),gj{s)), 

Jo Jo 

rT f-T 

^1,3=1 {F{un,h{^)),gj{s))ds, In,j= {an{un,h{^))h{s), gj{s))ds. 
Jo Jo 

Since /' G L^{[0,T]) and for every X G L'^{Q), {t,Lo) ^ ipjX{Lo) f'{t) G L'^i^T^H), 

(i) above implies that as n — > oo, /^ • — > J^ {uh{s),ipj)f'{s)ds weakly in L^(r2). Similarly, 

(v) implies that as n — > oo, /^ • — > J^ {Fh{s),gj{s))ds weakly in L'^{Q), while (vii) implies 

that I^j -^ /g [Sh{s),gj{s))ds weakly in L3(0). 

To prove the convergence of /^ •, as in [34] (see also [14]), let Vt denote the class of 

predictable processes in L^ (Qt , Lq{Hq, H)) with the inner product 

{G,J)vj,='k[ {G{s),J{s)) ds = E [ trace{G{s)QJ{s)*)ds. 
Jo '^ Jo 

The map T : Vt ^ L^{^) defined by T{G){t) = f^ {G{s)dW{s),gj{s)) is linear and 
continuous because of the Ito isometry. Furthermore, (vi) shows that for every G G Vt, 
as n^ oo, (o-„(n„,h)n„,G)^^ -^ {Sh,G)-pT weakly in L'^{Q). 

Finally, as n — > oo, P^C = '"n/i(0) -^ C ^^ H and by (iv), {un,h{T),gj{T)) converges to 
{uh{T),gj{T)) weakly in L?'{^). Therefore, as n — > oo, (4.14) leads to 

{uh{T),^3)f{T)= {^,^3)+ [ {uh{s),^j)f'{s)ds+ [ {Sh{s)dW{s),gj{s)) 

Jo Jo 

+ [ {Fh{s),gj{s))ds+ [ {Sh{s),gj{s))ds. (4.15) 

on 



^0 



For 5 > 0, A: > 1, t G [0,T], let A G H\-6,T + 6) be such that ||/fc||oo = 1, /fc = 1 
{—6, t — ^) and fk = on (i, T + 5) . Then /^ -^ ^(-S,t) i^i L^, ^i^d /(- ~^ ~^t in the sense 
of distributions. Hence as A; — > 00, (4.15) written with / := fk yields 

0= {i-Uh{t),^j) + / {Sh{s)dW{s),^j) + / {Fh{s),^j)ds+ I {Sh{s),^j)ds 
Jo Jo Jo 

for almost all t G [0, T]. This relation makes it possible to suppose (after some modifi- 
cation) that Uh{t) is weakly continuous in H for almost all u; G fi. Now note that j is 



26 



I. CHUESHOV AND A. MILLET 



arbitrary and E J^ |5/j(s)|? ds < oo; we deduce that for < t < T, 

Mt) =^+ [ Sh{s)dW{s) + / Fh{s)ds + / Sh{s)ds G H. 
Jo Jo Jo 

Moreover Jq Fh{s)ds G H. Let / = 1(^_s^t+5)i using again (4.15) we obtain 



(4.16) 



Uh{T)=i+ Sh{s)dW{s)+ Fh{s)ds+ Sh{s)ds. 
Jo Jo Jo 

This equation and (4.16) yield that Uh{T) = Uh{T) a.s. 
Step 3: In (4.16) we stiU have to prove that ds ® dP a.s. on Jly, one has 
Sh{s) = a{uh{s)), Fh{s) = F{uh{s)) and Sh{s) = a{uh{s)) h{s). 
To estabhsh these relations we use the same idea as in [34]. Let 

V e X = L^{nT,n) n L^{n, L°^{[o,T], H)) nL^{nT,v) . 

Suppose that L2 < 2 and let < 77 < —3-^; for every t G [0, T], set 



r{t) 



2Ri + 2C^ \\vis)\\ti + Li + 2JL^\his)\o + —\h{s)\l 
■- V 



ds, 



(4.17) 



where Cj^ is a function of 1] such that (4.1) holds. Then almost surely, < r(t) < cxd for 
all t G [0,T]. Moreover, we also have that 

r G L\n,L°°{0;T)), e"^ G L°°(Ot), r' G L^^t), r'e"^ G L°°(0, L^((0,r)). (4.18) 

Weak convergence in (iv) and the property P„^ — > ^ in i/ imply that 



E(|^,(r)|2 e-(^)) - E|CP < liminf [E(|n„,,(r)|2 g-^^^)) - E|P„C|' 



(4.19) 



We now apply Ito's formula to |ii(t)pe ^^^' for u = Uh and u = Un,h- This gives the 
relation 

E(|u(r)|2e-'^(^))-E|u(0)|2 = E / e-'-^'U{\u{s)\^]-^ [ r'{s)e-''^'^\u{s)\'^ds, 

Jo Jo 

which can be justified due to (4.18) and the property |up G L^(ri,L°°((0, T)). Using 
(4.16), (4.2) and letting u = v + (u — v) after simplification, from (4.19) we obtain 

E / e-''(^) [ - r'{s){\uh{s) - v{s)f + 2{uj,{s) - v{s) , vis))} + 2{Fi,{s) , Uh{s)) 



+ \Sh{s)\l^ + 2{Sh{s) , Uhis))] ds < liminf X„, 



(4.20) 



where 



T 



Xn=E I e-'-(^) [ - r'{s){\un,h{s) - v{s)\' + 2{un,h{s) - v{s) , v{s)) } 
Jo 

+ 2{F{un^h{s)),Un,h{s)) + |o-„(n„,/,(s))n„||^ + 2(a-(n.„^,,(s))/i(s), Un,h{s))] ds. 
The inequalities in (4.1), (C2), (C3), and also (4.17) and Schwarz's inequality imply that 

-'■^'^[-r'{s)\un,h{s)-v{s)\^ 



Yr,. 



E / e" 
Jo 



+ 2{F{Un^h{s)) - F{v{s)), Un^h{s) - v{s)) + |cr„(u„,;i(s)) n„ - an{v{s)) n„| 

+ 2({a„(n„,fe(s)) -^^(^(s))} h{s),Un,h{s) - v{s)) ds < 0. 



(4.21) 



Furthermore, X^ = Yn + X]i=i -^ni '^i*^ 
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z, 



E / e 
/o 



-r(s) 



2r'(s)K,;,(s)) - 7;(s), ^(s)) + 2(F(n„,;,(s)), r;(s)) 

+ 2{F{v{s)),Un,h{s)) - 2(F(w(s)), i'(s)) + 2{an{Un,h{sWn , £^(^(s))l^ 

+ 2(a„(n„,;,(s))/i(s),r;(s)) +2(a(t;(s)) h{s),Un,h{s)) -2{Pna{v{s))h{s),v{s)) 



ds, 



E / e 
/o 



-r(s) 



2(o-„(u„,/i(s))n„, ([cr(u(s))n„ - a{v{s))) 



\Pna{v{s))nn\h 



ds. 



The weak convergence properties (i)-(vii) imply that, as n — > oo, Z^ — > Z^ where 



E 



T 



-^(^) [ - 2r'is) {uhis) - vis), vis)) + 2{Fhis),vis)) + 2{Fivis)),Uhis)) 



- 2{Fivis)),v{s)) + 2iShis) , ^iv{s))). + 2(5,,(s), i;(s)) 
+ 2(a(«(s)) /i(s), nft(s)) - 2(a(«(s))/i(s), ^(s)) ] ds. (4.22) 

As for Z^ we note that the Lebesgue dominated convergence theorem imphes that 

E /" e~''^'^\aivis))iUn-IdHo)\lods^O as n^oo. 
Jo 

Using once more the dominated Lebesgue convergence theorem, we deduce that 

rT 



Z'^^-El e~'^^'Maivis))\l^ds as n ^ oo. 
/o ^ 



(4.23) 



Thus, (4.20)-(4.23) imply that for any v e X, 

Ey^e-^WJ - r'{s)\uhis) - vis)\^ + 2{Fh{s) - F{vis)),Uhis) - vis)) 

+ \Shis) - aivis))\l^ + 2(^Shis) - aivis))his) , Uhis) - vis)^ }ds < 0. 

Let V = Uh £ X; we conclude that Shis) = (T(n/j(s)), ds (d dF a.e. For A G 
L°°iQT,V), set vx = Uh — Xv ; then it is clear that vx G X. Applying (4.24) to v := vx 
and neglecting \aiuhis)) - crivxis))\l^, yields 

rT 



(4.24) 

., v e 



E / e 
/o 



-rx{s) 



\V,is)\vis)\' + 2X\{Fhis) - Fivxis)),vis)) 



+ 



Shis) - d-ivxis))his), vis)j I 



ds < 0, 



(4.25) 



where rxis) is given by (4.17) with vx instead of v. Using (C3) we obtain 
E / e-'^^'^\i[aivxis)) - aiuhis))]his) , vis))\ds 







T 



< |A| E / \his)\o \iis)\ ( \jLMs)\ + V^2 Ms) 
as A ^ 0. Hence, by the dominated convergence theorem. 



ds^O 



limE 







^^^^UShis) - aivxis))his) , vis))ds 



E / e 
/o 



-ro(s) 



5/1(5) - aiuhis))his) , ^(s)!^^. 



28 I. CHUESHOV AND A. MILLET 

Furthermore, (4.1) yields for A / and s G [0,T] 

\{F{v,{s)) - F{u,{s)),v{s))\ < C \X\ [\v{s)f + \\v{s)f + \v{s)\' ||n,(.)||y . 
Thus we deduce as A ^ 0, 

E / e-'-^'-'HFf,{s) - F{vx{s)),v{s))ds ^E [ e-'-^'^^HF^is) - F{uh{s)),v{s))ds. 
Jo Jo 

Thus, dividing (4.25) by A > (resp. A < 0) and letting A — > we obtain that for every 

V G L'^{Qj', V), which is a dense subset of L'^{Qj', V), 



T 



[Fhis) - F{uh{s)) , v{s)) + {Shis) - aiuh{s))h{s) , i{s)) 



ds = 0. 



Jo 
Hence a.e. for t G [0, T], (4.16) can be rewritten as 

Mt) = i+ I <y{uh{s))dW{s) + / [F{uh{s)) + d{uh{s))h{s)\ ds. (4.26) 

Jo Jo 

Furthermore, (i)-(iii) imply that 

Ml \\uh{t)fdt) < supE / \\un,h{t)fdt<C{l + E\^\^), 

^Jo ^ n Jo 

e( sup \uhit)\*) < supEf sup \un,hit)\^)<C{l + E\C\*), 

^0<t<T n ^0<t<T ' 



Ie( / \\^h{t)\\^dt\ < supE /" \\ur,,h{t)\\ttdt<C{l + E\C\^). (4.27) 

^Jo ^ n Jo 

This completes the proof of (2.31). 

Step 4: Now we prove that Uh G C{[0,T],H) almost surely. We first note that (4.26) 
yields that e~^^Uh. G C{[0,T],H) a.s. for any 6 > 0. Indeed, since for 5 > the operator 
e"''^ maps HtoV and V to H, we deduce that e"^^ /^ F{uh{s)) ds belongs to C([0, T],H). 
Condition (C3) implies that e~ f^ cr(n/j(s)) h{s) ds also belongs to C([0, T],H). Finally, 
condition (C2) implies E f^ \e-^'^a{uh){s)\l ds < +oo. Thus f^e-^'^a{uh{s))dW{s) 
belongs to C{[0,T],H) a.s. (see e.g. [12], Theorem 4.12). Therefore it is sufficient to 
prove that 

limE J sup \uh{t) - e-^^Uh{t)\'^ \ = 0. (4.28) 

5^0 [o<t<r J 

Let Gs = Id — e~ and apply Ito's formula to IGsUhit)]"^. This yields 

ft ft 



\GsUh{t)\' = \Gse-2j\\GsUhis)\\'ds + 2I{t) + J\Gsaiuhis))\i^ds 

+2 / {B{uhis)) + Riuhis)) + aiuhis))h{s), Gjuhis)) ds, (4.29) 

where I{t) = J^ [Gs'7{uh{s))dW{s),GsUh{s)). By the Burkholder-Davies-Gundy and 
Schwarz inequalities we have 

/ rT ^ 1/2 

E sup |/(t)| < CE{ \GsUh{s)\^\Gsa{uhis))\l ds 

0<t<T \Jo J 



< \e sup |G5n,,(t)|2 + ^E / \Gia{uh{s))\lds. 

^ 0<t<T ^ Jo 
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Hence for some constant C, (4.29) yields 

E sup \GsUh{t)f<2\GsC\^ + CEf \Gsa{uh{s))\lds 

0<t<T Jo 
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+ 4E 



[Biuhis)) + Riuhis)) + aiuhis))his), GJuhis) 



ds. 



Since for every u & H, \Gsu\ ^ as 5 — > and sup^^g \G5\l{h,h) ^ 1) '^s deduce that if 
{(fk} denotes an orthonormal basis in H, then \Gso'{uh{s))Q^''^ipk\^ — > for every k and 
almost every {uj,t) G i7 x [0,r]. Since 

sup\Gsa{uh)\l^ < J^sup|G5(TK)Qi/Vfc|' < C\a{uh)\l^ G L\n x [0,r]), 



<5>0 



5>0 



the Lebesgue dominated convergence theorem implies that E J^ |G5cr(u/i(s))|? ds — > 0. 
Furthermore, given u £ V we have ||G|u|| ^ as 5 — > and sup^^g IG^I^j^y^^) < 2. Hence 
(^B{uh{s)) + R{uh{s)) + a{uh{s))h{s), G'juh{s)) -^ for almost every {io,s). Therefore, as 
above, the Lebesgue dominated convergence theorem concludes the proof of (4.28). 

Step 5: To complete the proof of Theorem 2.4, we show that the solution u^ to (4.26) is 
unique in X := C([0, T],H)n L^{[0, T], V). Let w E X be another solution to (4.26) and 

TAT = inf{t > : \uh{t)\ > N} A mi{t > : \v{t)\ > N} AT. 

Since |u/i(-)l s-'^d \v{.)\ are a.s. bounded on [0, T], we have tn ^^ T a.s. as A^ ^ oo. 
Let U = Uh — V. By Ito's formula we have 

tAr pt^^N 

^-af, ''\\nHir)\\i,drp^^^^^^^2 ^ / q,{s)ds + <^{t A Tn), (4.30) 







where 



*(s) 



-o-Io \\^h{r)\\J^dr 



a\\uHis)fn\Uis)\' 



- 2\\U{s)f - 2{B{u,{s)) - B{v{s)), U{s)) + \a{n,{s)) - a{v{s))\l^ 
+ 2{[a{uh{s)) - a{v{s))]h{.s), U{s)) - 2{R{uh{s)) - R{v{s)) , U{s)) 



and 



$(r) = 2 / e 
/o 



-^Io\\M-)\\n'ir(^U{s),[a{uh{s))-a{v{s))]dW{s)). 



Now we set a = 2Crj where Cr^ is defined by (2.8). Then using (2.8) and Conditions (C2) 
and (C3) we obtain that for some non negative constant G{t]) which depends on -q, Ri, 
Li, Li, i = 1,2, and is independent of L2, 



+ (2R^ + L^ + ^\his)\l + 2 VLi|/i(s)|o )\Uis) 



< e 



-a Jo \\-'J-h{r)\\ii^dr 



(2 - 3r? - L2)\\U{s)f + C7(7?) (l + |/.(s)|g) |C/(s)p . (4.31) 



First consider the case of a general (random) control function h. Below we use the notations 



X{t) = sup 

0<s<t 



|e-«/o^"^ \\^hi'')\\ndr\u{s A r7v)p} , Y{t) 



tArjsi 



^Io\WH(r)\\^n'ir\\U(^S^2^^_ 
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Then it follows from (4.30) and (4.31) that for 3r/ < 2 - L2, 

X{t) + (2 - 3r? - L2)Y{t) < C{r]) / (l + \h{s)\l) X{s)ds + I{t), 

Jo 

where I(t) = supo<s<t |'5(s A TAr)|. An argument similar to that used to prove (4.10), 
based on the Burkholder-Davies-Gundy inequality, (C2) and Schwarz's inequality, yields 
that for t G [0, T] and /3 > 0, 

1/2 



E/(t) < 6E 



tArjv 

e 



-2a J^ ll"^MII«'^-|^(s)| V(n/.(s)) - a{v{s))\lds 



<PEX{t) + ^ / EXis)ds + ^EY{t) 
P Jo P 



Now we are in position to apply Lemma 4.1. If we choose t] = 1/3, 2(3 = exp{— C(l/3)(^+ 
M)}, then (4.3) holds under the condition L2 (1 + 36exp{2C(l/3)(r + M)}) < 1. There- 
fore, since supo<s<T | e""^o ^ II"''('')IIh'^''|C/(s A r7v)P [ < 2N, relation (4.4) implies that 
'EX{t) = for all t and hence, 

E sup |e-"/o^"'^ll«heOli«'^'^|C/(sAr;v)P|=0. (4.32) 

Since limAr^oo ''"TV = T a.s., and by (4.27) we have a.s. Jq ll'"/i('5)||[>^c?'5 < 00, we deduce 
that \U{s,Lu)\ = a.s. on 17^. Thus, we conclude that Uh{t) = v{t), a.s., for every t E [0, T] 
which yields the uniqueness statement in Theorem 2.4 for a general control function. 

Suppose now that we only have L2 < 2 and that h possesses a deterministic bound 
V'(t) G L^(0,r); let r] e]0, ^^]. Then it follows from (4.30) and (4.31) that 

VN{t) < C{ri) \l + |^(s)|2j VNis)ds with V^(t) = Ee^'^i'o^'^ \\u^{r)\\i,drp(^^ ^ ^^)|2_ 

Jo L -I 

Since the function s 1— > \ip\'^ belongs to L^{0,T), we can apply the Gronwall lemma to 
obtain (4.32) and to conclude the proof for the case considered. 

Finally suppose that a = a where a satisfies condition (C2) with L2 < 2. For h ^ 
set Wl^ = ^i + /o H^) ds and let P be the probability defined on {0,,J^t) by 

^=exp[-l\{s)dWs-\l\h{s)\Us). 

The Girsanov theorem implies that W is a IP Brownian motion with the same covariance 
operator Q. The above arguments prove that under P, the evolution equation 

Mt)=^+ f F{uh{s))ds+ [ a{uh{s))dW^ 
Jo Jo 

has a unique solution in X. Thus the Girsanov theorem implies that under P, (2.29) has 
a unique solution in X. Finally, once well-posedness is proved, computations similar to 
that used to obtain (4.5) in the case p = 2 can be used to deduce that (2.31) holds; this 
completes the proof of Theorem 2.4. □ 

Note that it follows from the consideration above that we only need the requirement 
(2.27) concerning the growth of a in order to obtain weak compactness of the sequence 
^n{'Un,h{s))h{s) in L'p{Qt]H) for p = 4/3 > 1. This weak compactness makes it possible 

to pass to the limit in the term E /^ {^n{un,h{s))h{s), v{s))ds in the expression for Z^ for 
elements v from the class X which contains the limiting function Uh- 
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